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Abstract 

The influence of motion of ions and electron temperature on nonlinear 
one-dimensional plasma waves with velocity close to the speed of light in 
vacuum is investigated. It is shown that although the wavebreaking field 
weakly depends on mass of ions, the nonlinear relativistic wavelength essen- 
tially changes. The nonlinearity leads to the increase of the strong plasma 
wavelength, while the motion of ions leads to the decrease of the wavelength. 
Both hydrodynamic approach and kinetic one, based on Vlasov-Poisson equa- 
tions, are used to investigate the relativistic strong plasma waves in a warm 
plasma. The existence of relativistic solitons in a thermal plasma is predicted. 

PACS numbers: 52.35.Mw, 52.35.Fp 
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I. INTRODUCTION 



Strong plasma waves passing through a plasma with phase velocity slightly 
smaller than the velocity of light are the subject of great interest during last 
two decades. Such waves can be excited in plasma by relativistic bunches of 
charged particles or laser pulses. The excited plasma waves can be used both 
to accelerate charged particles and to focus charged bunches [1]. Plasma- 
based accelerator concepts are currently under intensive development (see 
overview in Ref. [2] and numerous references therein). Accelerating gradi- 
ents in the plasma wave can reach the values of tens GeV/m (notice that in 
conventional linacs accelerating gradients are in order of tens MeV/m), that 
is confirmed in recent experiments [2]. Also the focusing fields can be much 
more than that reached in conventional focusing magnetic systems. The ac- 
celeration of charged particles by relativistic strong waves also is considering 
as a possible mechanism of ultra-high energy (up to 10 20 eV) cosmic ray 
generation in astrophysical plasma. 

In a cold plasma an amplitude of one-dimensional plasma wave is lim- 
ited by the wavebreaking field. In nonrelativistic case, when the wave phase 
velocity v p h is much less than velocity of light (v p h <C c), the wavebreaking 
amplitude is equal to [3] E* = m e cu pe v p h/\e\, where tu pe = (47rn e 2 /m e ) 1 / 2 
is the electron plasma frequency, n is the density of electrons in unper- 
turbed plasma, m e and e are the electron rest mass and their charge; the 
ions assumed to be immobile. In relativistic case the wavebreaking field is 
equal to [4] E rel = [2( 7 - l)} 1/2 /(3, here (3 = v ph /c, 7 = (1 - p 2 )~ 1 ' 2 is the 
relativistic factor and E re i is normalized on E*. The one-dimensional rela- 
tivistic strong waves (RSW) can be excited in a plasma by wide relativistic 
bunches of charged particles or intensive laser pulses [2] (when k p a ^> 1, 
where k p = u pe /v p h, a is the characteristic transverse sizes of bunches or 
pulses). 

Another important characteristic of nonlinear plasma waves is the depen- 
dence of the wavelength on the wave amplitude. Both in linear case and 
in nonlinear nonrelativistic one, the plasma wavelength in cold plasma is 
X p = 2nv p h/ oo pe . In the relativistic nonlinear regime, when plasma electrons 
get a relativistic velocity in the process of oscillations, nonlinear wavelength 
increases with the amplitude [5-8]. In the ultra-relativistic case (7 3> 1), for 
the wave amplitude E mp = E reh the wavelength approximately is equal to [5] 
4(2 7 ) 1 / 2 A p . One can see that for large 7 the wavelength is essentially more 



2 



than usual linear plasma wavelength. 

In the previous studies the plasma ions, in the process of oscillations, in- 
cluding the nonlinear relativistic waves, usually were assumed to be immobile 
due to their large mass. In Ref. [9] it is shown that when 7 <C (M /lQm e ) 1 ^ 
(here M is mass of an ion; for example, for hydrogen plasma, consisting of 
protons and electrons, this condition gives 7 < 5), the wavebreaking ampli- 
tude approximately is equal to E re i and the motion of ions can be neglected. 
However, the dispersion properties of the relativistic strong plasma waves, 
which take into consideration the motion of ions, are not elucidated up to 
now. This problem has been considered in Sec. II on the base of cold hy- 
drodynamics equations for the arbitrary 7 and mass of particles forming the 
plasma. The necessity to take into account the ion motion conditioned by 
the following reasons. Firstly, because the maximum relativistic wavelength 
and amplitude grow proportional to 7 1 / 2 , the plasma ions (even heavy ions) 
in such strong field can get velocity which is sufficient to give essential con- 
tribution in the process of charge separation in the wave. On the other hand, 
in a semiconductor plasma positively-charged particles (holes) have a mass 
similar or less than that of electron. The problem has also astrophysical 
aspect. The pole region of the pulsars is considered to be filled with an 
electron-positron plasma in which the strong plasma waves can be excited 
and high energy charged particles generated [10]. It is obvious that in the 
plasma wave passing through an electron-positron plasma, neither electrons 
nor positrons may be considered as neutralizing background. 

Another important problem is the influence of plasma temperature on 
RSW. Finite plasma temperature has decisive significance for description of 
RSW near "wavebreaking". Actually, according to one-dimensional theory 
of relativistic plasma waves in a cold plasma, at the "breaking" point hy- 
drodynamic plasma electron velocity is equal to the phase velocity and the 
density of electrons n e tends to infinity [5-7, 11] (in this case spatial behavior 
of the density is similar to the ^-function). On the other hand, in a thermal 
plasma (even when the temperature is low), the pressure tends to infinity 
when n e — > 00. Thus, in this case, the pressure, as well as the plasma tem- 
perature, should be taken into consideration. In Ref. [12] the finite plasma 
temperature effect on nonrelativistic (v p h <C c) wavebreaking field is consid- 
ered using ID waterbag model for the distribution function. In this model 
it is assumed, that the electron distribution function during oscillations is 
constant in a limited interval of velocities and is equal to zero outside of 
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this interval. It is shown [12] that maximum amplitude of the plasma waves 
decreases with the temperature. In Ref. [13] the ID relativistic waterbag 
model is used to investigate RSW in a warm plasma. Using relativistic equa- 
tion of motion with the pressure term for plasma electrons, in Ref. [11] it 
is shown that in the case 7 3> m e v ph /3T (where T is the temperature of 
electrons) the wavebreaking field is proportional to T -1 / 4 . The authors of 
Ref. [14] have analyzed the influence of low temperature (T <C m e c 2 ) on 
excitation of nonlinear wake fields by a relativistic charged bunches. They 
considered the equation of motion obtained using second moments of the 
distribution function. In the present paper (Sec. Ill) the hydrodynamics 
equations are used to study dispersion properties of RSW in a warm plasma. 
The dispersion correlation for weakly nonlinear case is obtained. In Sec. IV 
the strong plasma waves are investigated on the base of relativistic Vlasov 
kinetic equation and the Poisson equation. 

II. ION MOTION EFFECT ON DISPERSION PROPERTIES 
OF RELATIVISTIC STRONG PLASMA WAVES 

In this section we consider a cold uniform plasma consisting of positively- 
charged particles (for example, protons or positrons) with mass m + and 
electric charge q + , and negatively-charged particles (electrons or negatively- 
charged ions) with mass m_ and the charge q_. The relativistic equation 
of motion and the continuity equation for each plasma component and the 
Poisson equation for one-dimensional steady plasma waves are: 

dz \q_\ 

p^-W±M = 0, (2) 
dz dz 

^ = l-N„ + \q + /q„\N + , (3) 

where z = k p (Z — v p ht), k p = co p /v ph , uj p = (47rn -^/m_) 1 / 2 , n _ is the 
density of negatively- charged particles in equilibrium, (3± = v±/c are dimen- 
sionless velocities, 7± = (1 — Z^) -1 / 2 , densities N± are normalized on the 
equilibrium values. The electric field strength is normalized on the nonrela- 
tivistic wavebreaking field m-0J p v p h/\q-\ and obeys the formula 
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E(z) = -(l/(3 2 )d<S>/dz, (4) 

where $ = $_ = 1 + \q_\(p/m_c 2 > I/7, (p is the electric potential. From 
expressions (1), (2) and (4) we have: 

(3 ± = [(3-^l-T 2 ) 1,2 ]l{P 2 + ®l\ (5) 

iV ± = M$ ± /($|-7~ 2 ) 1/2 -/?]. (6) 

Substituting N±($ ± ) and expression (4) in the Poisson equation (3) one 
obtain the following differential equation of the second order for $: 

d ^ + BW\ *± ? 1 = (7) 

Here $ + = f — q + Lp/m + c 2 = 1 + — $) > I/7 and /x = |g + /g_|m_/m + . 
The electric potential f is assumed to be equal to zero when plasma density 
is equal to the equilibrium density. 

Equation (7) can be rewritten in the form 

d 2 <S> dU 

+ ^ = 0. (8) 



dz 2 d$ 

U = /3 3 7 2 {[/3-(^ + -7- 2 ) 1/2 ]/^ + [/3-($ 2 -7- 2 ) 1/2 ]}- 

Here, for convenience, U(<&) is chosen to be equal to zero in a point $ = 1, 
where it reaches a minimum. When fx — > 0, Eq. (7) reduces to the known 
equation for nonlinear waves in a plasma with immobile ions [2,7]. Formally, 
Eq. (8) describes one-dimensional motion of a particle in a field with potential 
[/(<&); the values $ and E correspond to the coordinate and velocity of this 
fictitious particle respectively. Function U determines the characteristic of 
the field moving through the plasma. In Fig. 1 this function is presented 
for 7 = 10 (f3 ~ 0.995) and different values of fx. One can see that for the 
arbitrary parameters the solutions of Eq. (8) [or Eq. (7)] are the periodic 
plasma waves (including the wave with zero amplitude-unperturbed plasma). 
Integrating Eq. (8) we have 

^ = -p 2 E = ±[2(U nmx -U)] 1 / 2 , (9) 
az 
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where C/ max is maximum value of U(&) in the process of oscillations. From (9) 
it follows that the plasma wave amplitude is equal to E mp = (2f/ max ) 1 / 2 / / 9 2 . 
Substituting maximum permissible value of U(<&), which reaching at $ = 
I/7, in this expression, we find the wavebreaking field: 

E WB = 2 1 / 2 7 [l + (l-e i 1/2 d /2 )//i], (10) 
6 = 6 = l + /i(7-l)/(7 + l). 

In the case /^C 1, from (10) follows the expression 

^^(l + ^/8)[2( 7 -l)] 1/2 /A (11) 

which reduces to the well known relativistic wavebreaking field for plasma 
with immobile ions, when (j, — [4]. Fig. 2 shows the wavebreaking field E WB 
depending on ji (for example, for electron-positron plasma \i = m e /m pos = 1, 
for the hydrogen plasma /i = m e /m prot ps 5.455 x 10~ 4 ) for different values of 
7. Both in nonrelativistic case and in relativistic one the wavebreaking field 
weakly increases with \x. For example, according to (10), in nonrelativistic 
case (7 ~ 1), E WB (/j, = 1) only 2(1 — 2 -1 / 2 ) 1 / 2 pz 1.08 times exceeds the 
wavebreaking amplitude at fx — 0. Proceeding from the shape of the "poten- 
tial" [/(<&) (see Fig. 1) one can expect that the plasma wavelength undergoes 
considerable change with fi. In Fig. 3 the dependence of relativistic plasma 
wavelength A p [note, that according to the variables accepted in (l)-(3), the 
linear plasma wave at /i — corresponds to the value A p = 2n] on amplitude 
presented. The curve 1 corresponds to the case of immobile ions and coincides 
with that previously obtained [7]; in this case A p grows with the amplitude 
due to relativistic velocities of the oscillating plasma particles (notice, that 
in nonlinear nonrelativistic regime the plasma wavelength does not depend 
on amplitude). The motion of positive ions for fixed amplitude causes the 
decrease of charge separation length and therefore, leads to the decrease of 
the wavelength. Fig. 3 clearly shows competition of two tendencies. For the 
small \x the wavelength grows with the amplitude due to nonlinearity. With 
the increase of E mp , the effect of ion motion becomes more essential. When 
fi is not small, the behavior of the wavelength caused mainly by ion motion. 
In electron-positron plasma (/i = 1), A p monotonously decreases with the 
increase of E mp , in contrast to the case of heavy ions (// w 0). The results 
of simulations presented in Fig. 3 conform with the well known result of 
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linear theory (E mp <C 1; see, e.g., Ref. [15]): A p = A p0 /(1 + f^) 1 ^ 2 , where 
A p o = A p (/i = 0). Note also that the results practically did not change for 
arbitrary 7 ^> 1. Considerable decrease of the relativistic plasma wavelength 
with \i is demonstrated in Fig. 4. 

Previous studies have shown that energy of electrons (or positrons) ac- 
celerated in the field of relativistic nonlinear wave, in a cold plasma with 
immobile ions, can reach a value of 4m e c 2 7 3 [7,16]. In the general case the 
relativistic factor of a resonant electron passing from a point with the dimen- 
sionless potential $1 to a point with $ 2 is equal to [7] 

7ac C ~7a CC (0)+2 7 2 ($ 2 -$ 1 ), (12) 

where 7 acc (0) ~ 7 is value of the relativistic factor at the initial point. When 
pi — 0, <3> min = I/7 < $ < $ max ps 27 [7]. Substituting this maximum and 
minimum values in (12), for maximum energy of accelerated electrons one 
obtains (7 acc )max ~ ^l 3 - With \i growth, the maximum energy decreases due 
to the decrease of $ max (see Fig. 1). For the case of electron-positron plasma 
function U ($) is symmetric with reference to axis $ = 1 and, as it is easy to 
see, in this case $ m i n = 1/7, $ ma x = 2 — I/7. Then, the maximum energy of 
accelerated particles is (7 aC c)max ~ 47 2 , that is in 7 times less than that in 
the case /i — 0. 

If a bunch of charged particles with density n^z) and electric charge q b 
passes through a plasma, adding to the left side of equation (7) the value 
a(z) = P 2 (qb/\q-\)n b (z)/n _, we obtain the equation that describes the ex- 
citation of steady plasma wake fields by the bunch. In this case the phase 
velocity is equal to velocity of the bunch. In this section the properties 
of RSW have been investigated by simulation of wake wave generation by 
charged bunches. 

Above we have considered the case // < 1. However, one can see that the 
obtained results are valid also for fi > 1, if E we replace by — E [this new 
E is normalized to m + uj p v p h/q + , tu p = (47m 0+ q , +/m + ) 1 / 2 ], jx replace by l//i, 
replace the subscript " + " by " — " and vice versa. 

III. THE INFLUENCE OF ELECTRON TEMPERATURE ON 
RELATIVISTIC NONLINEAR PLASMA WAVES: 
HYDRODYNAMIC APPROACH 

Here we continue to consider the dispersion properties of RSW in the 
frame of hydrodynamic approach, investigate relativistic nonlinear waves in a 
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warm plasma. Adding the relativistic pressure term — ('y 2 /N e )(l — (3(3 e )dP/dz 
[11,14] with P = r(A jr e /7 e ) 3 [11] (which is relativistic generalization of usual 
equation of state for one-dimensional adiabatic compression) to the equation 
of motion of plasma electrons one can obtain the equations: 



= + <13) 

dE 1 d 2 § 

1-N e . 



dz (3 2 dz 2 



In (13) f3 e = v e /c and r = T/m e c 2 are dimensionless velocity and temperature 
of plasma electrons, 7 e = (1 — /J 2 ) -1 / 2 , $ = 1 + \e\ip/m e c 2 . The plasma ions 
are assumed to be immobile due to their large mass. The density of electrons 
N e = n e /n normalized to the unperturbed value no, as usually, is obtained 
from the continuity equation: 

N e = P/(P-p e ). (14) 

When t — > equations (13) and (14) describe RSW in a cold plasma [5,6]. 
Note, that the value r = 1 corresponds to the temperature of about 6xl0 9 K. 
For laboratory plasmas the temperature changes in the bounds r ~ 1CT 6 -j- 
10~ 2 ; for star plasmas r ~ 10~ 5 — 1. 

The dispersion correlation can be obtained analytically for weakly non- 
linear wave, when u = f3 e (z)//3 <^ 1. In this case from equations (13) and 
(14) we have: 



d 2 u ( du 

(a Q - a x u + a 2 u 2 )-^ - (a x - 2a 2 u) \ — \ + u + u 2 + = 0, (15) 



a = l-3r/(3 2 , ai = l + 3r(3-2/3 2 )//3 2 , a 2 = 3/3 2 /2-3r(6-ll/3 2 /2 + /3 4 )//3 2 . 
Looking for solution of Eq. (15) as (see, e.g., Ref. [17]) 

u = eux(*) +e 2 u 2 {^) +e 3 u 3 {^) + 
d^/dz = \p/A p = k + eki + e 2 n 2 + 
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we obtain 



A p /A p = a 1/2 (l + b(3 2 J, (16) 
b = _3/16 + 3ra 1 (10-9/3 2 + /5 4 )/8/3 4 + 3r 2 a 2 (2-/5 2 ) 2 //9 6 , 

where e = f3 m / f3 <C 1 is the small parameter, (5 m = (/3 e ) max and A p is 
the wavelength. In the linear case {(3^ — > 0) from (16) it follows that the 
wavelength decreases with the temperature according to well known Bohm- 
Gross dispersion correlation: A p = A p (l — 3T//3 2 ) 1 / 2 . On the other hand, 
in a cold plasma (r = 0) the wavelength increases due to nonlinearity (see, 
e.g., Ref. [15]): A p A p (l + 3/3^/16). Fig. 5 shows the dependence of the 
wavelength on wave amplitude in thermal plasma obtained by simulation of 
Eqs. (13) and (14). In the case of low temperature this dependence almost 
coincides with that in a cold plasma (compare curves 1 in Figs. 5 and 3). 

In the case (3 — > 1 equations (13) and (14) can be easily integrated (see 
also Ref. [11]): 



$ - 




0. (17) 



One can see that the thermal term can not be neglected near the "wavebreak- 
ing" (f3 e — > 1) even for low temperature. In the latter case the wavebreaking 
field is proportional to r^ 1 / 4 [11]. 

In the frame of hydrodynamic theory the velocity of plasma electrons can 
not exceed the wave phase velocity. Actually, if (3 e > (3, then, according 
to expression (14), the density of electrons becomes negative, that has no 
physical sense. In reality, when j3 e ~ (3 in a warm plasma (even when the 
temperature is low), considerable part of electrons gets velocities more than 
the phase velocity due to their thermal energy distribution. When the tem- 
perature is not low, the energy distribution effect on plasma waves is essential 
in all cases. Therefore, the strong waves near the "wavebreaking" and when 
r is not low, can be described correctly in the frame of kinetic approach. 

IV. KINETIC THEORY OF THE RELATIVISTIC STRONG 

PLASMA WAVES 

As in the previous section, here we assume plasma ions to be immobile. 
The kinetic approach, for one-dimensional steady fields passing through a 
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warm plasma, gives the following system, obtained from relativistic Vlasov 
equation and Maxwell equations (see, e.g., Ref. [18]): 



P- 



V 

(l+p2)l/2 
~dz* 



dz 



<9$ df 

dz dp 



0, 



+ (3 2 (1-N e ) = 0, 



+oo 



f(p,z)dp, 



(18) 
(19) 
(20) 



where p = p z is the plasma electron momentum, normalized to m e c, f(p, z) is 
the distribution function. In an unperturbed plasma the distribution function 
is equal to ID relativistic Maxwell distribution [19] 



fo = A[l + (l+p 2 Y/ 2 /T}eM-a+P 2 ) 1/2 M, (21) 
A = t/2K 2 (1/t), 
(p 2 ) = r[K 1 (l/r)/K 2 (l/r)+4r], 



((1 + p 2 ) 1/2 >o = 2r + (1 - r 2 )K 1 (l/r)/^ 2 (l/r), 

where K n (x) is the modified Bessel function of the n-th order. In (21) we 
also have written out average squared pulse and total energy for the one- 
dimensional equilibrium distribution / , that may be interesting for future 
investigations. 

Solving Eq. (18) by the method of characteristics (see, e.g., Ref. [20]) 
and requiring that the function / reduces to the equilibrium distribution (21) 
at $ = 1, one obtains the following general solution: 



/ = A(l + S/r)eM-SM, (22) 
S = (l + 9 2 ) 1/2 , S = -7 2 [/?r±(r 2 - 7 ~ 2 )n 
r = /5p-(l+j9 2 ) 1/2 + $-l, 

In the expression for g, the plus sign corresponds to the case p < and the 
minus sign to p > (3% in equilibrium ($ = 1) we have g — p. Substituting 
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the expressions (20) and (22) in Eq. (19), one obtains equation for $. The 
plasma electron density A^ e ($) obtained numerically from (20) and (22) is 
presented in Fig. 6. In the case of low temperature, the integral in (20) can 
be calculated by the Laplace asymptotic method [19]. The value of N e is at 
maximum when $ pa I/7 and is equal to 

iV max « [r(l/4)/4] 7 (/? 7 A) 1/2 (2/r) 1 / 4 » O.M^) 1 ^" 1 / 4 . (23) 
According to expression (23), in a cold plasma iV max — > 00, that conforms 
with the previous investigations [5,6,11]. When r C 1 and $ > I/7, the 
dependence iV e ($) approximately is described by expression (6). 

Simulations of the problem show that the plasma wavelength increases 
with the wave amplitude and tends to infinity for a solitary wave (soliton). 
Fig. 7 shows dependence of the maximum value of the amplitude (which 
corresponds to the maximum of electric field strength in the soliton E s ) on 
plasma temperature for different 7. One can see that E s is almost constant 
and equal to the relativistic wavebreaking field E rd for the values of r up to 
0.05 — 0.1 and then decreases rapidly. 

As it was mentioned above, the relativistic plasma waves can be excited 
by charged bunches or laser pulses. In order to describe the excitation of the 
wake field by a charged bunch, it is necessary to add in the left side of Eq. 
(19) the value (3 2 a(z) [the definition of ct(z) see in Sec. II]. The nonlinear 
periodical wave and the solitary wave excited by uniform electron bunch are 
plotted in Fig. 8. The plasma electron density behind the soliton tends 
to its equilibrium value (iV e — > 1) and the electric field strength tends to 
zero. However, in this case the average plasma electron momentum (p) = 
f-m pf(p, z)dp/ f(p, z)dp tends to a non-zero constant value. This does 
not seem strange because the solitary wave can be considered as a wave 
with infinite wavelength. Hence, bulk motion of plasma electrons behind 
the solitary wave takes place, while the plasma remains neutral. Equations 
(18)-(21) have also non-periodical solutions. However, such solutions have no 
physical sense [19] and should be considered in the frame of non-stationary 
kinetic theory. Thus, in a thermal plasma with immobile ions two kinds of 
steady waves can exist: periodical waves and solitons. 

V. CONCLUSIONS 

The results presented in this paper supplement the theory of nonlinear 
relativistic plasma waves, taking into consideration motion of ions and finite 
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plasma temperature. It is shown that the nonlinearity leads to the increase 
of relativistic wavelength, while ion motion leads to the wavelength decrease. 
For example, in electron-positron plasma the wavelength monotonously de- 
creases as the amplitude increases. The relativistic wavebreaking field weakly 
depends on the ion mass. 

Contrary to the case of cold plasma, in a warm plasma the relativistic soli- 
tary waves (solitons) can exist. The plasma wavelength grows monotonously 
with the amplitude in the warm plasma due to nonlinearity. It has been 
found that maximum electron density in the plasma wave decreases with the 
temperature as T" 1 / 4 and tends to infinity in a cold plasma, that was shown 
by previous investigations. 
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Figure 1. "Potential" U{<P) for the value 7=10. 
1- (i=0; 2- (i=0.1; 3- (x=1. 
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Figure 2. The wavebreaking field depending on ji. 
1-7=1.01; 2-y=1.5; 3-7=3; 4-7=10. 




Figure 3. Relativistic plasma wavelength A p depending on the electric 
field amplitude E mp ; y=1 0. 1 - (i=0; 2- (i=0.05; 3- (i=0.1 . 4- ^i=0.5; 5- (i=1 . 




Figure 4. Relativistic strong wave in plasma with immobile ions (1- (i=0) 
and in electron-positron plasma (2- (i=1); y=10. 




Figure 5. The relativistic plasma wavelength in a warm plasma as a 
function of the wave amplitude; y=1 0. 1 - x=0.001 ; 2- x=0.01 ; 3- x=0.04. 




Figure 6. Plasma electron density as function of the dimensionless electric 
potential 4> for y=3. 1-t=10" 4 ; 2-t=10" 3 ; 3-t=10" 2 . 
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Figure 7. The maximum value of electric field amplitude of nonlinear 
plasma wave depending on plasma electron temperature. 
1- 7=3; 2- 7=6; 3- 7=10. 
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Figure 8. Electric field strength (a), plasma electron density (b), and average 
electron momentum (c) in strong plasma field, excited by uniform electron bunch 
(7=3, t=0.1). The rectangles show the bunch. 
1- periodic wave, rit/n^OA; 2- solitary wave, ni/no=0.6575. 



